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Abstract 

We introduce and study algebraic dynamical systems generated 
by triangular systems of rational functions. We obtain several results 
^ ■ about the degree growth and linear independence of iterates as well 

as about possible lengths of trajectories generated by such dynamical 
systems over finite fields. Some of these results are generalisations of 
those known in the polynomial case, some are new even in this case. 
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1 Introduction 

Let F be an arbitrary field F and let F±, . . . , F m G F(Xl, . . . , X m ) be m 
rational functions in m variables over F. For each i = 1, . . . , m we define the 
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k-th iteration of the rational function Fj by the recurrence relation 

F{ 0) =X t , Ft^F^Ft^,-..^^), £ = 1,2,.... (1) 

In this paper we consider dynamical systems generated by multivariate ratio- 
nal functions, we refer to [1, 21, 22] for a background on algebraic dynamical 
systems. 

More precisely, we define the vectors u„ = (u ni i, . . . , M n ,m) G F m by the 
recurrence relation 

u n+ i,i = F i (u n> i,...,u n ,m), ra = 0,l,..., i = l,...,m, (2) 

with some initial vector Uq = (uo,i, • • • , uo, m ) G F m . 

As we work with rational functions, we make the standard convention 
(see [3, 9, 10]) that 

(T 1 = 0. (3) 
Using the following vector notation 

F = (F^Xi, . . . , X rn ), . . . , F m (X±, . . . , X m )), 

we have the recurrence relation 

u„+i=F(u„), n = 0,l,.... (4) 

In particular, for any n > and i — 1, . . . , m we have 

«n,i = ^(Uo) = Fj; n \uo,l, • • • , l*o,m) 

or 

u n = F^(uo), 

provided that u„ has been generated by (4) without using the convention (3) 
(that is, no poles have been encountered). 

Clearly, if we work over a finite field of q elements, the above sequence (4) 
of vectors {u n } is eventually periodic with some period r < q m . 

One of the important characteristics of the dynamical system generated 
by Fi, . . . , F m G F(Xl, . . . , X m ) is the degree growth of the functions (1). It 
is of great interest for the theory of dynamical systems and has been studied 
in a number of works, see, for example, [2, 25] and references therein. It is 
also important for applications to pseudorandom number generators [24]. 
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More precisely, although for a "typical" system an exponential degree 
growth is expected, there are several examples of systems where the degree 
grows much slower (which is highly beneficial for their applications), and 
such systems are of special interest. 

For example, in [14, 17] several types of multivariate polynomial systems 
T = {F 1 , ... , F m } of m polynomials in m variables over a finite field ¥ q have 
been constructed and studied, having the "triangular" form 

F 1 (X 1 , . . . , X m ) = X 1 G 1 (X 2 , . . . , X m ) + H 1 (X 2 , . . . , X rn ), 

F m -i(Xi, . . . , X m ) = X m -iG m -i(X m ) + H m -i(X m ), 
F m {Xii • • • ■> X m ) g m X m + h m , 

with d, Hi e ¥ q [X i+1 , . . . , X m ], i = 1, . . . , m - 1, and g m , h m e F„ g m ^ 0. 
These systems have been further investigated in [13, 18, 19, 20]. 

For the systems (5), in the case of constant polynomials G{ G F* in [14] 
and polynomials Gj with leading terms of special form in [17, 18], a series 
of results have been obtained about the distribution of the corresponding 
sequences given by (2) that are much stronger than those known for generic 
systems. Moreover, for these classes of polynomials, it has been shown in [17] 
that the degrees of the iterations of the polynomials Fj, % — 1, . . . , m, grow 
significantly slower than the exponential growth expected for the iterations 
of a "generic" system of m polynomials in m variables. In turn, this leads 
(see [18]) to much better estimates of exponential sums, and thus of discrep- 
ancy, for vectors generated by (5) than for those originated from arbitrary 
polynomial systems (see [4, 5, 16]). 

We also note that the results obtained in [17, 18] regarding the degree 
growth of the iterations of the polynomials in (5) hold over any field F. 

In this paper we extend the class of rational dynamical systems with slow 
degree growth and present an analogue of the construction (5), but with 
rational functions defined by 

Fi(Xi, . . . , X rn ) = Xi 1 Gi(X 2 , . . . , X m ) + Hi(X 2 , . . . , X m ), 

(Xl, . . . , X m ) — X e ™_^G m -i(X m ) + H m -i(X m ), 
F m (Xi, . . . , X m ) = g m X^ + h m , 

with e 1 ,...,e m e {-1, 1}, G h Hi e F[X i+1 , . . . , X m ], i = 1, . . . , m - 1, and 
g m , h m e F, g m ^ 0. 
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We note that for m—1 and e = 1 we obtain the classical linear congru- 
ential generator which have been successully used for decades in the theory 
of Quasi Monte Carlo methods, see [7, 8], and for m — 1 and e = —1, the 
classical inversive generator, see [9, 10, 11, 12]. 

For the above class of multivariate rational functions, we study the de- 
gree growth under iterations and, using an approach similar to that of [17, 
Lemma 1], we show in Section 3 that under certain additional conditions 
imposed on the systems of rational functions (6), the degree grows polyno- 
mially. 

Moreover, for applications to pseudorandom number generators, follow- 
ing the standard technique almost identical to that of [17], one almost im- 
mediately obtains bounds on the exponential sums with elements of the se- 
quence (4) generated by the system (6) (satisfying the conditions outlined in 
Section 3), that in turn leads to estimates on the uniformity of distribution 
of the vectors (4). However, one has also to prove that for any k ^ I and 
nonzero vector a = (ai, . . . , a m _i) G F m_1 , the linear combination 

m—1 

Q^E^-if) (7) 

1=1 

is a non-constant rational function. We note that in the case of rational 
functions this does not follow directly from the degree argument as in the 
case of the polynomial systems (5) in [17], but we give such a result in 
Section 4. 

Since the derivation of such bounds of exponential sums for our systems 
does not bring anything new to the area, we do not do this here but rather 
concentrate on the study of the degree and linear independence of iterates, 
which is also of interest for the general area of algebraic dynamics. 

Furthermore, we consider a related question about the length of trajecto- 
ries generated by iterations (4) over a finite field ¥ g . We remark that in this 
case a trajectory falls into a cycle if u t = u s for some integers t > s > 0. In 
particular, we show that under some rather broad conditions for any fixed 
e > 0, for all but o(q m ) initial vectors u G F™, the trajectory length t of the 
iterations (4) is at least q 1 ^^ 6 . 

We note that Silverman [23] has considered a question about periods of 
general polynomial systems but in somewhat dual situation when the initial 
value is fixed and the iterations are considered over a family of finite fields. 
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The results of [23] apply to very general systems, however the estimates are 
only logarithmic rather than a power of the field size. 

Moreover, we give necessary and sufficient conditions for the systems (6) 
to generate sequences of maximal period. We note that for the case = 1 
for alH = 1, . . . , m, the maximal period length of the sequence generated by 
the system (5) is achieved whenever the conditions of [15, Theorem 6] are 
satisfied. Our result is a generalisation of that of [15]. 

2 Structure of the Iterations 

As in [17], we can describe explicitly the iterations of the rational functions 
Fi as follows. 

Let us define the sets 

1+ = {1 < i < m : d = 1} and i_ = {1 < i < m : e { = -1}. (8) 

We also define 

G?\x t+1 ,...,X m ) = ^(F^,...,^- 1 )), 
Hf\x i+u ...,X m ) = H,{F^\...,Ft^. 

Lemma 1. Let Fi,...,F m be rational functions defined by (6). Then for 
i = 1, . . . , m — 1 and k = 0, 1, . . for the rational functions F^ given 

by (l), 

1. for every i G I + , i < m, we have 

F^ = XiG i)k + H ijk , (9) 
where G ijk , H i)k e ¥(X i+1 , . . . , X m ) are defined by 



^i,k — LrjtJo . . . Lr. 



(2) n {k) 

H hk = HiG^p ...Gf ] + H^Gf . . . Gf ] + ... + H^G^ + H?\ 

(10) 

2. for every i e I-, i < m, we have: 

F (k) = XiR itk + S^k 
XiRi tk -i + Si,fc-i 
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where Ri,k,Si,k are defined by the recurrence relations 

tti,k — -tt«,fc-2 + tl i tLi y k~l 

C o i zj( k ) Q 

>Ji,k — Lrj Ji,k-2 + n i Di t k-1 

for k > 1, with the initial rational functions 

Rifi — 1) Si,o — 0, Ri t i = Hi, S it i = Gf, 
3. ifm<El + , then 

= g k m X m + (g^ 1 + ... + g m + l)h rn ; 

4- if m e then 

F(k) = (A k ) hl X m + (A% 2 
m ~ (A") 2tl X m + (A") 2 y 

h m g m \ k = ( {A k ) ltl (A k ) 1>2 
1 J - \ (A k ) 2<1 (A k )2,2 



(12) 



where 

A k = 



Proof. The case = 1, i = 1, . . . , m, is given by [18, Lemma 1]. We con- 
sider now that = — 1 and prove the result by induction on the number of 
iterations k. For k — 1 it is clear from the definition of the system, so we 
consider the statement true for the first k — 1 iterations and we prove it for 
the A;-th iteration. For i — 1, . . . ,m — 1, we have 

F (k) _ F ( F (k-i) F (k-i) F (k-ih 

F (k-i) H (k) G (k) XR k _ 1+ s Hi k) (fc) 

pi(k—l) ■Xi-Rj,fc-i+5'i,fc-i 

and thus we conclude this case. When e m = — 1, it is also clear as the /c-th 
iteration of 

_ h m X m -\- g m 
fm — 77 
A m 

is given by A k as simple calculations show. □ 



6 



We want to describe the degree growth of the iterations of the rational 
functions defined by (6), and in particular to prove that we have the same 
effect of slow degree growth as for the polynomial systems (5) described in [17, 
Lemma 1]. To be able to give an explicit formula for the degree growth we 
need to impose some further conditions on the degrees of the polynomials G{ 
and Hi, i — 1, . . . , m — 1. 

Let Fi,...,F m be rational functions defined by (6). From now on we 
consider the system (6) satisfying the following conditions for Fi for any 
i — 1, . . . , m: 

1. if ei = 1, as in [17, 18], we assume that the polynomial Gi has a unique 
leading monomial X^^ 1 . . . Xm' m , that is 

where & G F* and Gi G ¥[X i+1 , . . . , X m ] with 

deg x . Gi < Si j, deg x . Hi < s id , j = i + 1, . . . , m; (13) 

2. if Ci = — 1, we assume that the polynomial Hi has a unique leading 
monomial X%\ +1 . . . X^ m , that is 

H i = h i X£? 1 ...X%<» + H i , 

where hi G F* and H { G W[X i+1 , . . . , X m ], and 

deg x . flj < Sij, deg x . G { < 2s id , j = i + 1, . . . , m. (14) 

We note that having these conditions also allows us to consider the ra- 
tional function system with constant multipliers Gj, i — 1, . . . ,m — 1. We 
remark that in [14], the case of constant polynomials Gi, i — 1, . . . ,m — 1, 
in the system (5) was considered, but this case is different from the case 
of rational functions as the conditions on the degrees also differ, see (13) 
and (14). Having this, we prove the following formula for the degree growth 
which coincides with [17, Lemma 1]. 
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3 Degree Growth 



Theorem 2. Let Fi,...,F m be rational functions defined by (6) satisfying 
the conditions (13) and (14) and such that s^+i ^ 0, z = 1, . . . , m — 1. T/ien 
i/ie degrees of the iterations of Fi, . . . , F m grow as follows 

degi^ (fe) = - — ^-r T A: m ^Si i i + i...s m _i im + ^(/i;), z = 0, . . . , ro - 1, 
(m — 1)\ 

deg^W = 1, 

where ifti(T) G Q[T] is a polynomial of degree degipi < m — i. 

Proof. The proof is based on Lemma 1. The case = 1, using (9) and (11), 
follows exactly the same as in [17, Lemma 1]. 

We prove now the case when e» = — 1. Using the conditions (13) and (14) 
and the recurrence relation (12), it is easy to see that 

deg F t {k) = deg Ri, k + 1 = deg H[ k) + 1 

= deg HiH® ...H^ + 1 = J2 deg + I. ^ 

3=1 

As in [17, Lemma 1], we use induction on the number of variables m. For 
m = 2 we easily see that degif^ = degiJi = si j2 , and thus, by (15), we 
have that degF^ = ks lt2 + 1. We assume now that the theorem is true for 
m — 1 variables and we prove it for m. For any i = 1, . . . , m — 1, by the 
induction hypothesis, we have 

k k 

deg Fj k) = ^g + 1 = £ deg /^(i^, . . . , F^) + 1 



As 



E d ^g ((Jfe 1 ')*" 1 ■ ■ ■ (^T 1 ^) + i 

3=1 
k ( 1 

/ , ( 7 : 7Tr(j _ I) 7 " 1 s i,i+lSi+l,i+2 ■ ■ ■ S m -l,m+ 

V (m — i — 1 ! 

...+(j-l)s + 1- 

fc 1 

Vr 1 " 1 = — ^^(/c + 1) - B m -i(o)), 



3=1 
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where B m _i is the Bernoulli polynomial of degree m — i (which has the leading 
coefficient equal to 1), we finally obtain the desired result. □ 

4 Linear Independence 

Theorem 3. Let F 1 , . . . , F m be rational functions defined by (6) satisfying 
the conditions (13) and (14) and such that s^+i ^ 0, i = 1, . . . ,m — 1. Then, 
for k 7^ I and a nonzero vector a e F m_1 ; Qk,i, a is a non- constant rational 
function. 

Proof. The proof reduces to proving that deg(Fi^ — fP) > 1, where s < 
m — 1 is the smallest index such that a s ^ 0, as the variable X s does not 
appear in the polynomial 

m— 1 

Qfc,/, a - a s (F^ - F^) = ^(if'-jf). 

i=s+l 

If e s = 1, it is clear, as from Theorem 2, for k > I we have deg G St k > deg G s j. 
If e s = — 1, by (11), we have 

p(k) _ p(l) _ X sRs,k + S^k X s R s j + S s> i _ Uk^s 

X s R S) k-i + S St k-i X s R St i_i + S Sj i-i Vk,i, s 

where 

Uk,i, s = Xg(R Sj kR s ,i-i — R s ,k-iRs,i) 

+ X s (R Si kS s> i-i + S s> kR s ,i-i — R s ,k-iS s ,i — S s> k-iR s ,i) 
+ S s ,kS s ,i-i — S Sj k-iS S: i. 

and 

Vk,i, s = XgR St k-iR s ,k-i + X s (R St k-iS s j-i + R St i-iS Si k-i) + S St k-iS s ,i-i- 

Without loss of generality we may assume that k > I. Using (12), we 
obtain 

R s ,kR s ,i-i — R s ,k-iR s ,i = (G^R s ,k-2 + H^R St k-i)R s ,i-i 

— R s ,k-i(G^ R s ,i-2 + H^Rgj-x), 
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and thus, using Lemma 1 and Theorem 2, we derive 
deg( J R s , fcJ R Si/ _ 1 - Rs^Rsj) = deg(H^ - # 

(16) 

= degtff 'Rs^Rs^ > degi2 a , fc _ii2 8) i_i > 1 
for k > I, which concludes the proof. □ 
Note that, as in [18], we can include m-term linear combinations 

m 

QM,a = £^ (fe) -*f) 
i=l 

with a G F m , but in the case of a\ — . . . — a m _i = 0, a rn ^ 0, the nontriviality 
also depends on the divisibility of k — I by the multiplicative order of g m . 



5 Trajectory Lengths 

In this section we work over a finite field ¥ q . 

Theorem 4. Let Fi,...,F m be rational functions defined by (6) satisfying 
the conditions (13) and (14) and such that s^+i ^ 0, i = 1, . . . ,m — 1. Then, 
for any T > 1 for all but 0{T^q m ~ 1 ) initial vectors u G F™, the trajectory 
length of the iterations (4) exceeds T. 

Proof. Let U be the set of u = . . . , u m ) e such that 

Gi(iii+i, ...,u m ) = 

for some z = 1, . . . , m — 1. Clearly #W = 0(g m_1 ). 

We now build a sequence of sets £4, A; = 0, 1, . . ., recursively. 
We put U = U. 

Assume that Uq, . . . ,Uk are defined and let 

w fe = w u...uw fc . 

Then we let U k+ i be the set of the initial values 

u e \ m 

such that for the corresponding sequence of vectors (4) we have u fc+1 e U. 
Inspecting (6), we now see that, by our assumption, for any v G F™, there 
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is a unique preimage u G F™ \ W under the map given by (6) (that is, with 
v = F(u)). In turn, we see that for any v£W there is a unique corresponding 
initial value u G F™ \ V\4 with u^+i = v. Thus jftAk = jflA. 

Since there are O^q™" 1 ) vectors u G F™ that contain a zero in at least 
one component, we see that the set £t of initial values for which, for some 
integer t <T, the vector u t has a zero component, satisfies 

S T = 0{q m ~ x + T#W) = 0(Tq m - 1 ). (17) 

We see that if a vector u G F™ \ £ T generates a trajectory of lengths 
t <T then u t = u s for some nonnegative integer s < t. 

Now, if e m _i = 1, then we remove the set Tt of initial vectors u G F™ 
such that 

G m -i,t(uo) = G , m _i )S (u ) 

for some integers s and t with T>t>s>0. By Lemma 1 we see that 
G m -i,t — C m _i jS is a nontrivial polynomial of degree 0(t). Hence, 

#jr T = o( t( l m ~ l ) = {^q™- 1 ) ■ (18) 

\0<s<t<T J 

Furthermore, if e m -\ — — 1, then we remove the set J-r of initial vectors 
u G F™ such that 

^m-l,t( u o)'Rm-l,s-l( u o) — ^m-l,t-l( u o)-Rm-l,s( u o) 

for some integers s and t with T > t > s > 0. As in the proof of Theorem 3 
(in particular, see (16)) we note that Lemma 1 implies that R m -i,tRm-i,s-i — 

-i,t-iRm-i,s is a nontrivial polynomial of degree 0(t). Hence, again we 
obtain the bound (18). 

We remark that for T > t > s > 0, for any solution u = (wo,i, • • • , «o,m) G 
F™ \ T T to the equation 

fW(u )=fW(u ), 

the component Wo.m-i is uniquely defined by u 0jm . So there at most q m ~ x 
such solutions for every fixed t and s with T > t > s > and thus at 
most T 2 ^"^ 1 for such t and s. Combining this bound with (17) and (18) we 
conclude the proof. □ 
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Clearly if the map u h-> F(u) is a permutation, as for example, in [13], 
then all trajectories are purely periodic. So we always have s = in the argu- 
ment of the proof of Theorem 4. This leads to a better estimate 0(T 2 q m ~ v ) 
on the number of initial values generating trajectories of length at most T. 

6 Maximal Periods 

In this section we show that the periods of the rational function systems 
over F q defined by (6) with e$ = — 1 for all i = 1, . . . , to are given by the 
orbit lengths of certain linear fractional transformations, also called Mobius 
transformations. In particular, we describe the case when the systems (6) 
achieve maximal periods in their orbits. We also note that in [15, Theorem 
6] there are given necessary and sufficient conditions for the system (6) to 
achieve maximal period in the case = 1 for all i — 1, . . . , m. 
We denote 

Uo,i = («o,i+i, ■ • • , «o,m) e F™"\ i = 1, • • • , m - 1. 

Lemma 5. Let F ± , . . . , F m e F 9 [Xi, . . . , X m ] be as in (6) with = —1 /or a// 
i = 1, . . . , to. Assume that the sequence generated by the lower m — i rational 
functions F i+1 , . . . , F m in F™~* zs purely periodic with period r i+i for some 
i = 1, . . . , to — 1 . Then we have the following description for the kri + ±-th 
iteration of Fi on any initial vector u G F™: 

F^ +l) (u ) = /f ) K,), k>l, (19) 

where Ri, n+l and Si tTi+1 are defined by (12) and fi is the Mobius transforma- 
tion in the variable Y , 

( u o,i) ' 

In particular, the orbit length of Fi in u is given by the orbit length of fi in 

U ,i- 

Proof. We first note that the orbit length of Fi is a multiple of Tj+i. Indeed, 
let Ti be the orbit length of the system F i: . . . , F m in the initial vector u . 
Then Tj = lcm (r i+ i, 7/j), where r\i is the period ofthe sequence {u nti } defined 
by the iterations of the polynomial F iy and thus is a multiple of r i+1 . This 
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shows that, in order to describe the period of F, . . . , F m on the initial vector 
u , it is enough to consider only kr i+ i-th iterations of Fj. 
By (11) we have 

F (Tl+l) (uo) = u ^ iRi ' Ti +^ n °^ + S i,n+i( u o,i) = j.( Mq .) (20) 

Mo,i-Ri,T i+ i-l(Uo,i) + <Si,T i+ i-l(Uo,i) 

Now, 

F/^ +l) (uo) 

=F (n + >) (Ff-^Vo), ^r 1)Ti+l) (u ), • • • , il (fe - 1)T ' +l) (uo)) 

^^(Ff-^^uo),^,...,^) (21) 
^((fc 1)T * +l) (uo)-Ri, Ti+1 (uo,i) + -S i)Ti+1 (u ,i) 



F ((k 1 ) Ti + 1 )( Uo ) j r. iTj+i _ 1 ( Uo ,.) _|_ S'j^j-^uci) 

To prove (19) we use induction over k. For k — 1 it is clear. 

We now assume that the statement is true for k — 1 and we prove it for 
k. Using (20), (21) and the induction hypothesis, we derive 

* ( o) = 7m7 — ^ r~ ^ ; c = ji { 

Jj l«0,iJ-Ki,r i+ i-HU 0l iJ + <3 i)7 . i+1 _i^U ,i; 

which concludes the proof. □ 

Lemma 6. Let T = {F ± , . . . , F m } be a system of polynomials over¥ q defined 
by (6). Let i — 1, . . . , m such that = — 1 in the system (6). Then we have 

-H-i,k>Ji,k-l — -fti,fc-l<->i,fc — i — J-J ^i^rj ■ • • (Jj , 

where Ri,k,Si,k are defined by (12). 
Proof. We use induction on fc. If A; = 1, by (12), we get 

Ri,i£>i,o ~ RifiSi^i = —G{. 

We assume the statement true for k and we prove it for k + 1. By (12) and 
the induction hypothesis we have 

Ri,k+lSi,k ~ Ri,kSi,k+l 

—(G\ k+1 ^ Ri,k-i + H\ k+1 ^ Ri : k)Si,k — Ri,k{Gf +1 ^ Si^i + H- k+1 ^S iik ) 
— — g\ + \Ri tk Si^-i — Ri,k-\Si,k) 
=(-l) k+1 G i G<?\..G\ k) G ( i k+1) 
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and thus we conclude the proof. 



□ 



As usual, we say that a polynomial / G ¥ q [X] of degree d > 1 is primitive 
if it is the minimal polynomial over ¥ q of a primitive element of ¥ g d (that is, 
an element of multiplicative order q d — 1), see [6]. 

Next, we present necessary and sufficient conditions for the system (6) to 
achieve maximal period over the prime field F p . 

Using [15, Lemma 2] which holds for the functions Fi in the system (6) 
for which = 1, we have the following analogue of [15, Lemma 5] (with an 
almost identical proof which we do not present here). 

Lemma 7. Let T = {F±, . . . , F m } be a system of polynomials over ¥ p defined 
by (6). Let the index 1 < i < m such that e, = 1 and assume that the period 
of the sequence generated by the lower m — i polynomials F i+1 , . . . , F m in ¥ p n ~ t 
is p m ~ % and that G i:P m-i(u 0ji ) = 1. Then, for the rational functions H i:pm -i 
defined by (10), we have 

H ijP m-i(u 0ti ) = ^2 ^i( V )> 

where 

R l = H t G?K.. G r i) . (22) 

Now, using Lemma 7 and [15, Theorem 6] for Fj with e { = 1 in the 
system (6), we have the following result. 

We recall that the sets I + and I are given by (8). 

Theorem 8. Let T = {F 1 ,...,F m } be a system of polynomials over ¥ p 
defined by (6). Then the sequence {u n } generated by (4) is purely periodic 
with period r = p m if and only if the following conditions are satisfied 

1. for every i G I + , i < m, we have 

II Gi{v) = l and ^( v )^°; 

2. for every i G i_ ; i < m, we have: 

(a) if Ri :P m-i _i{\Xq) = 0, then 

R ip m-i{\\Q) = S iiP m-i _ 1 (uq) and S i:P m-i(uQ)S ij pm-i_i( y u ) 7^ 0; 
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(b) if Ri^m-i -i{u Q ) ^ 0, then 



X 2 _ Rj^prn-ijUo) _ IlveF™- 1 G l( V ) 



i?j p m-i_ 1 (Uo) i2j jp m-i_ 1 (Uo) 

is a primitive polynomial over ¥ p ; 
3. if m G i+, £/ien g m = 1; 

^. ifm&I-, then X 2 — /i m X — g m is a primitive polynomial over F p . 

Proof. We prove the result by induction on m. For m = 1, if m G i+, then 
it's clear that the period p is achieved if and only if g m — 1 and /i m G F*. If 
m 6 L, we have 

Fi = ^Xf 1 + h u 

which, by [3, Theorem 1], has maximal period p if and only if the polynomial 
X 2 — h\X — gi is a primitive polynomial over F p . 

We assume the statement true for the first m — 1 variables and we want 
to prove it for m. Let the sequence {u n l } = {(u nj 2, ■ ■ ■ ,\ m )} be defined 
by the last m — 1 components of the vectors in the sequence {u„}. By the 
induction hypothesis we know that the period r of the sequence {u n ,i} is 
p m_1 , and taking into account the first remark in the proof of Lemma 5, we 
see that the period of {u n l } is of the form kp™" 1 , for some 1 < k < q. 
Thus, proving the maximality of the period of {u„} reduces to proving that 

(k) 

k = q. We note that the representation of F x given by Lemma 1 does not 
depend how we choose e2, . . . , e m in the functions F 2 , . . . , F m , but only the 
functions Gi t k,Hi,k or G^\h^ if 1 e /+ or 1 e respectively. Thus, the 
representation of given by Lemma 1 is the same regardless if ij G I + or 
ij G /_ for 2 < j < m. 

Thus, the case 1 G I + follows identically as in the proof of [15, Theorem], 
and we do not repeat it here. 

We consider now the case 1 G i_. By Lemma 5 we have 

F^ m - 1 \u ) = ff(u , l ), k>l, 

where _ _ 

j. — ^F liP m-i(u 0il ) + *Si, p ™-i(uo,i) 

yi? 1)P m-i_ 1 (u 0i i) + 5 , i, p ™-i-i(u 0i i) 
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and thus the maximal period of the sequence generated by the iterations of 
Fi is given by the case when fi achieves maximal orbit length in w ,i- 

We distinguish now two cases. If i? 1)J ,m-i_ 1 (u 0) i) = 0, then we note that 
S'i ) pm-i_i(uo J i) 7^ 0, as otherwise fi = 0. Thus, we have the case of linear 
generator 

j — ^^ZljL^ll y 1 l ^ 1 .p m " 1 ( u o,i) 



S , i )P m-i_ 1 (uo,i) <S'i,p m - 1 -i(uo,i) 
The maximal period is achieved in this case if and only if 

■R lip m-l(u 0[1 ) _ ^ S 1 , p m-l(UQ A ) _^ ^ 

S'l ) pm-l_ 1 (Uo ! l) 5'l ) pm-l_ 1 (Uo ! l) 

which is equivalent to 

i2 ljp m-l(Uo,l) = 5' 1) pm-l_ 1 (Uo i l), lS'l j pm-l(Uo,l) ^ 0. 

We consider now the case of i?i ) pm-i_ 1 (u 0)1 ) 7^ 0. We note that in this 
case fi achieves maximal period if and only if, under a linear transformation, 
it has the same property. Taking into account that i?i,pm-i_ 1 (u ,i) 7^ 0, we 
can make the linear transformation 

Y — > i?l ; pm-l_ 1 (Uo,l) 1 Y — i?l jp m-l_l(U0,l) 1 »S'l i pm-l_l(Uo i l) 

and obtain the following inversive generator which, by a slightly abuse of 
notation, we denote also by }\, 

n /^x — -Ri,p'«- 1 (uo,i)'S'i i pm-i_i(uo,i) + ^l,p"i-l-l(u 0l l)5'l ) pm-l(Uo i l) x 

^ Y > ~ r — ^?r^ 



i? lp m-l_ 1 (u 0i l) 

Applying now Lemma 6 we have 

i?l ) pm-l(Uo,l)S'l ) pm-l_ 1 (Uo,l)— -Rl,pm-l_ 1 (Uo,l)5'l j pm-l(Uo,l) 

= -G 1 (u OA )G?\u OA )...G ( f~ 1 \uo,i)- 

Let T = {F 2 , . . . , F m }. Now, as the period induced by T is p m_1 , the 
elements 

UOA.^UOA),---,^" 1 ^^!) 
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are all the distinct elements of F™ *, and thus we obtain 



v 



G 1 (u , 1 )GY>(u , l )...G? >(u ^) 

= G 1 (G 0il )G 1 (^(ii 0il ))...G 1 (^ 1 - 1 )(ii 0il ))= J] G ^)- 



This concludes that 



= n vg F;- lg i( v ) rl + Ri, pm - i(u o ,0 



i? 1)P m-i(u 0i i) -Ri, P m - 1 -i(uo,i) 

Applying now [3, Theorem 1], we know that fi achieves maximal period q if 
and only if the polynomial 

yigti) ru^ gi(v) 

— t~ — 7 — t~ — r 

i? lp m-l_ 1 (u 0j l) i?l j pm-l_ 1 (u 0) l) 

is a primitive polynomial over F p . □ 
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